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Abstract 

We study the average density of resonances (DOR) for a semi-infinite disordered 
chain, coupled to the outside world by a (semi-infinite) perfect lead. A set of equa- 
tions is derived, which provides the general framework for calculating the average 
DOR, for an arbitrary disorder and coupling strength. These general equations are 
applied to the case of weak coupling and an asymptotically exact expression for the 
averaged DOR is derived, in the limit of small resonance width. This expression is 
universal, in the sense that it holds for any degree of disorder and everywhere in 
the (unperturbed) energy band. 



1 Introduction 

Open quantum systems often exhibit the phenomenon of resonances. Resonances corre- 
spond to quasi-stationary states which have a long life-time but eventually decay into the 
continuum. (A particle, initially within the system, escapes to infinity.) One approach 
to the problem of resonances is based on the study of the analytic properties of the scat- 
tering matrix, S(E), in the complex energy plane. Resonances correspond to the poles, 
E a = E a — |r a , of S(E) on the unphysical sheet [1, 2]. An alternative approach amounts 
to solving the stationary Schrodinger equation, with the boundary condition of an out- 
going wave only. This condition, which makes the problem non-Hermitian, describes a 
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particle ejected from the system. The Schrodinger equation with such boundary condition 
admits complex eigenvalues E a , which correspond to the resonances [1, 2]. This kind of 
approach leads in a natural way to a non-Hermitian effective Hamiltonian, and it has been 
used for a long time in scattering theory, including scattering in chaotic and disordered 
systems ([3, 4], and references therein). There are many examples of resonances in atomic 
and nuclear physics. More recently, there has been much interest in resonant phenomena 
in the field of chaotic and disordered systems (for a recent review see [4]). 

There is considerable amount of work concerning the distribution P(T) of resonance 
widths in one-dimensional disordered chains [5, 6, 7, 8, 9]. Numerical studies presented 
in that work demonstrate that, in a broad range of T, P(T) oc T~ 7 , with the exponent 
7 being close to 1. (A similar behavior pertains also to two- and three-dimensional 
systems with localized states [7, 9]). An analytical calculation [8] has been performed for 
a one-dimensional continuous (white-noise) potential, in the semiclassical limit, when the 
localization length is large. It has been shown in that work that in a broad range of T, 
P(T) is well fitted by a function r -1-25 . For sufficiently small T, however, (1/T) -behavior 
was obtained, followed by a sharp cutoff at still smaller T, due to the finite size of the 
sample. A different analytical approach was recently developed in [10], for a discrete 
tight-binding random chain. The limit of strong disorder (i.e., opposite to that of Ref. 
[8]) was considered and the (1/T) -behavior (for a semi-infinite chain) was derived. 

In the present paper we develop a new approach to the problem of resonances. The 
approach is based on counting the number of poles of a resolvent of the corresponding 
non-Hermitian Hamiltonian. We derive a set of equations for a semi-infinite disordered 
chain, coupled to a (semi-infinite) lead. This set contains, in principal, the full solution 
of the problem, for an arbitrary coupling strength and disorder. The equations simplify 
considerably in the weak coupling regime. For this case we rigorously derive the asymp- 
totically exact (l/r)-rule and present a simple scaling formula, which contains only the 
product of the localization length and the density of states, at the relevant energy. 

2 The model and its Effective Hamiltonian 

The system is depicted in Fig.l. 
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Fig. 1: A one-dimensional chain coupled to a lead. Sites of the disordered chain are denoted by 
black dots. Sites of the lead are denoted by open circles. The arrow represents the out-going wave. 
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Black dots, labelled by n — 1, 2, . . ., designate sites along the semi-infinite disordered 
chain. Each site of the chain is assigned a site energy, e n . Different e n 's (n — 1, 2, . . .) are 
independent random variables chosen from some distribution g(e). Open circles, labeled 
by n = 0, — 1, —2, . . . represent a perfect semi- infinite lead to which the chain is coupled. 
All sites of the lead are assigned e n = 0. The lead simulates the free space outside the 
chain. All nearest neighbor sites of the chain are coupled to each other by a hopping 
amplitude t, and the same is true for all nearest neighbor sites of the lead. The only 
exception to this rule is the pair (0,1) which provides coupling between the chain and the 
lead. The hopping amplitude for this pair is taken to be equal t'. This allows us to tune 
the coupling from t' — (closed chain) to t' — t (fully coupled chain). 

As was mentioned in the Introduction, the most direct approach to the problem of 
resonances is based on solving the stationary Schrodinger equation with the boundary 
condition of an outgoing wave (see Fig.l). The Schrodinger equation for the entire system 
(chain +lead) is a set of coupled equations: 



-tyn+i ~ tipn-i = Eip n (for n < 0) (1) 

-*V>_i - tVi = EiJ in = 0) (2) 

-^ 2 -tVo + ei^i = Efa (ra = l) (3) 

-tlp n+1 - t^ n -i + € n i) n = Elp n (n > 1) (4) 



We recall that e n = for n < 1 (the lead) and it is random for n > 1 (the chain). 
Eqs.(l-4) are to be solved subjected to the boundary condition of a outgoing wave in the 
lead, i.e. ip n oc exp(— ikn) , for n < 0, with Imk < . The complex wave vector k is 
related to E by E = —It cos k. Using the plane wave shape of the outgoing wave in the 
lead, it is straightforward to eliminate from Eqs.(l-4) all ipn's with n < 1 , thus reducing 
the problem to a system of equations for the amplitudes ip n on the sites of the disordered 
chain alone (n — 1,2,...): 

-tyn+x - ttpn^ + l n i) n = Ei\) n (n = 1, 2, . . .) (5) 

with the condition ip = 0. Here e n = e n for n — 2, 3, . . ., but not for n — 1. This end site 
is assigned a complex energy 

ei = ei — trjexpik, (6) 

where the parameter rj = (t'/t) 2 describes the coupling strength to the outside world. 
Thus, the resonances are given by the complex eigenvalues of the non-Hermitian effective 
Hamiltonian defined in (5). This non-Hermitian Hamiltonian, H, differs from the Hermi- 
tian Hamiltonian, H, of the corresponding closed system (i.e., with rj = 0) only by the 
complex correction to the energy of the first site (Eq.(6)): 

H — H — tr]e fk P, (7) 
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where P is the projection on site 1. We set the energy scale by taking t — 1 and denote 
by z the complex variable E/t. Note that, since k is related to z via z = —2 cos A;, the 
effective Hamiltonian H depends on z. Thus, Eq. (5) does not constitute a standard 
eigenvalue problem and the eigenvalues of H have to be determined self-consistently. 
Defining formally the resolvent G(z) = (z — H)' 1 we can write 

G — (z - H + rje^P)' 1 = [(z - H)(l + ^GP)}' 1 = 1 G, (8) 

1 + r\e lK GP 

where G = (z — H)^ 1 is the "unperturbed" resolvent. Resonances correspond to the 
singularities of the matrix G nm (z), in the complex z-plane, or to the roots of the equation 

l+ V e fk ^Gu(z) = 0, (9) 
where Gu(z) is the (1, l)-element of G in the site representation. 

3 An expression for the average density of resonances 

Let us write Gu(z) as 

g»w = ,- ei - gl( , r < 10 > 

where Si(z) is the self-energy for site 1. Eq.(9) can be written as 

F{z) = z-e 1 -S 1 (z)+r ] e fk{z) =0. (11) 

Denoting the zeroes of F(z) by z a , we write the density of resonances (DOR), for a given 
realization of the disorder, as p(z) = J2 a $( z — z a)- This expression refers directly to the 
semi-infinite chain and it should be understood as the N — > oo limit of the analogous 
expression for a finite chain, of size N. The limit is well defined for any Imz ^ and no 
division of the sum by iV is necessary [10], in contrast to the usual case of the density of 
states (on the real axis) for a Hermitian problem. 

It is convenient to define a ^-function of a complex variable F as 5(F) = 5(ReF)8(ImF) 
and to use the relation 

J25(z-z a ) = \F'(z)\ 2 8(F(z)), (12) 

a 

where the prime indicates a derivative with respect to z. This relation is a generalization 
of the corresponding relation for a real variable. It is valid if F(z) is analytic in a domain 
containing the zeroes z a and F'(z a ) ^ 0. Indeed, under such conditions in the neighbor- 
hood of a zero, the equation w = F(z) can be solved by an analytic function z = z(w) 
and the jacobian of the change of variables from (Rew, Imw) to (Rez, Imz) is |F'(^)| 2 
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due to the Cauchy-Riemann equations. We will assume that F'(z a ) 7^ with probability 
one. Substituting Eq.(ll) for F(z), we have 



p(z) = |1 - S[(z) + V X'(z)\ 2 5(z - ei - S^z) + V X(z)), (13) 

where X(z) stands for the expression e lk( ^ z \ 

Next we average Eq.(13) over all realizations of the set {e n } of the random site energies 
(we denote this average by angular brackets). Since the self-energy Si(z) does not depend 
on the energy e\ of the first site, we can single out the variable e\ and average over it 
explicitly. This leads to the following expression for the average DOR: 

<p ( z ) >= J deq{t) < [\l + r)X'\ 2 -2Re({l + r)X')S' 1 ) + \S[\ 2 }5{z-e-S 1 +r ] X) >, (14) 

The two correlated random variables, Si and S[ depend on the set {€2,63,...} of the 
random site energies. Note that we do not need to know the joint probability distribution 
for Si and S[. Indeed, only three special combinations of these two variables appear in 
Eq.(14). It is therefore useful to introduce three functions 

f(w) = <6(S 1 -w)> (15) 

g(w) = <S[(z)5(S l -w)> (16) 

h(w) = <\S[(z)\ 2 8(S 1 -w)> (17) 

and write (14) as 

<p(z)> = J deqie^l + rjX'iz^fiz-e + rjXiz)) 

-2Re(l + r]X'(z))g(z - e + r]X(z)) + h(z - e + r]X(z))}. (18) 

We emphasize that Si refers to a closed (semi-infinite) chain, so that one can use the 
standard relations for various quantities of the Hermitian problem [11]. In particular, 
we'll need the relation 

where 62(2) is the self energy for site 2, with site 1 excluded. It immediately follows that 

s^H-s^xi-s^)). (20) 

With the help of the definitions (15)-(17) and the recursion relations (19), (20) it is 
straightforward to derive integral relations between the functions /, g and h. For instance, 
using (20), h(w) can be rewritten as 

h(w) =< (1 - 2ReS' 2 + \S' 2 \ 2 )5{z - e 2 - S 2 - -) > . (21) 

w 
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The variable €2 is now "isolated" (i.e., all other quantities in (21) do not depend on this 
variable) and can be averaged over. This leads to 

h(w) = I deq(e) \f(z - e - -) - 2Reg(z - e - -) + h(z - e - -)]. (22) 
J w w w 

The two other equations can be derived in a similar way, yielding: 

9H = ^f deq(e)[g(z - e - -) - f(z - e - -)} (23) 
|w| 4 J WW 

and 

/H = nS / deq(e)f(z -e--). (24) 
|wp J w 

Let us stress that the above defined functions of the complex argument w are not complex 
valued functions in the usual sense but rather shorthand notations for a function of two 
real variables, w\ and w-i- For instance, Eq. (24) for the function f{w), which is the 
probability distribution for the real and imaginary part of the self-energy S± — w± + IW2, 
can be written more explicitly as 

f(w 1 ,w 2 ) = , 2] 2\2 1 de q( e )f( x - e — 2 Wl 2 ' v+ 2^ 2 2 )- ( 25 ) 

The three equations, (22)-(24), supplemented by the expression (18) provide the gen- 
eral framework for computing the average DOR, < p(z) >, in the complex plane z = x+iy, 
for any strength of disorder and for an arbitrary coupling 77. We do not attempt to solve 
the problem in its full generality but rather restrict ourselves to the small 77, i.e. weak 
coupling case. 



4 The weak coupling limit 

From now on we assume that the coupling constant rj is small and develop a "linear 
response theory" with respect to 77. In this limit the width of all resonances becomes 
proportional to r], so that the "cloud" of resonances in the (x,y)-plane gets squeezed 
towards the real axes. In order to define a meaningful i] — > limit for the DOR, one must 
stretch the y-axis by a factor of I/77. It is natural to switch from the original variable 
z = x + iy to a new variable, Z = X + iY, where X = x and Y = —y/r], where the minus 
sign accounts for the fact that the resonances are located in the lower half plane of the 
complex variable z. Thus, in the Z-plane resonances appear in the upper half-plane and 
Y is the resonance width, in units of rj. The limiting (i.e., in the 77 — > limit) average 
DOR, < p(X,Y) >, in the complex Z-plane is related to the original DOR, < p(x,y) >, 
as: 

< p(X, Y) >= lim 77 < p(X, -rjY) > . (26) 
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Next, we turn our attention to the functions /, g and h in the 77 — > limit. The 
function f(wi,w 2 ) is the probability distribution for the real and imaginary part of the 
self-energy Si. The shape of this function depends on the point z = x + iy. In the 77 — > 
limit, as has been just explained, one should set y = —r]Y, where Y does not depend on 
77. Furthermore, since the imaginary part of Si is proportional to 77, we set w 2 = rjW2- 
The distribution, F(W\,W2), for the new variables {W\ = wi,W2 = i]~ 1 W2) is related 
to the distribution f(wi,W2) for the "old" variables by F(Wi,W2) = 77/(Wi, 77W2). The 
distribution F(Wi, W 2 ) has a well defined r\ — > limit and it satisfies an integral equation, 
obtained from (25) by transforming to the new variables and taking the 77 — > limit: 

= wt I deqie)Fix -'-k' - Y+ wt ] - (27) 

The simplifying feature of the 77 — ^ limit is that functions g and h are easily ex- 
pressible in terms of /. More precisely, we have to define functions G and H, of the new 
variables: 

G(W 1 ,W 2 )=r l g(Wi,r l W 2 ), H(Wi, W 2 ) = vKW u V W 2 ). (28) 

These functions, in complete analogy with F(Wi,W 2 ), have a well defined 77 — > limit 
and satisfy integral relations which are derived from (22), (23) by writing them in explicit 
form (compare to (25)), transforming to the new variables and taking the i] — > limit. 
This yields: 



-Y + tttI) 



and 



r 1 W9 

H(W U W 2 ) = j deq(e)[H(X-e-—, -Y + ^) 



-2^G(X- e -^, -y + ^) + F(X-e-^, (30) 
One can easily check that (29) and (30) are satisfied by 

G(W 1 , W 2 ) = Y-F(Wi, W 2 ), H{Wi, W 2 ) = (^) 2 F(W U W 2 ) (31) 

Now we can derive an expression for < p(X,Y) > (p(X,y)}, in terms of the function F, 
by taking the 77 — > limit in Eq.(18) and using (26) and (31): 

< p(X, Y) >= ( /m ^ X ) )2 J rfeg ( e ) F (x - € ,-Y + ImX(X)). (32) 
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Note that, since z = X + ir]Y, in the 77 — limit k(z) approaches k(X) which is related 
to X by X = —2 cos k. Therefore Im\ in (32) is given by 

(33) 

It is convenient to define a new variable, J = W2/Y, so that J is the imaginary part 
of self-energy (the resonance width) in units of Yrj. Finally, we denote W\ by R (the 
real part of the self-energy) and define the probability distribution P(R,J), instead of 
F(Wi, W2). P(R, J) satisfies the integral equation 

P(R, J ) = jrJ dR '<i( x ~ \ ~ R 'W, ^ - i), (34) 

where the integration variable e was replaced by R' = x — e — R~ x . The expression (32) 
for the average DOR now reads: 



<p(X,Y)>=^(l-^)j dRq(X-R)P(R, I^l-^-l). (35) 

While performing integration in equations (34), (35), one should keep in mind that the 
function P(R, J) is identically zero for R > \fj . This property follows from the basic 
recursion relation for the self energy S. 

The integral equation (34), supplemented by the expression (35), completely defines 
the problem of resonances in the weak coupling limit. Note that the coupling strength, 77, 
does not appear explicitly in these equations. It only determines the units in which the 
resonance width is measured (the widths of all resonances is proportional to 77). Thus, 
< p(X, Y) > is determined solely by the properties of the closed system and, in this 
sense, Eqs.(34), (35) describe the regime of linear response with respect to the coupling 
strength rj. The function P(R, J) describes the joint probability distribution for the real 
and imaginary part of the self-energy (of the closed system and for z approaches the real 
axis) and the integral equation (34) has appeared previously in the study of Anderson 
localization [13]. Integration of (34) over J yields the integral equation for the probability 
distribution, P(R), of the real part of the self-energy: 

P{R) = ±j dR'q(x - 1 - R')P(R'). (36) 

This equation is very useful in the theory of one-dimensional localization [11, 12], because 
the knowledge of P(R) allows one to compute the localization length, £, according to: 

— 1 —= f dR'P(R')ln\R'\. (37) 

We shall use this relation in the next Section. 
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5 Large- J asymptotics 



We were not able to obtain the complete analytic solution of the integral equation (34). 
It is possible, however, to find the large- J asymptotics of the function P(R, J). We state 
the result and give the proof later: To leading order in the small quantity, 1/ J, 

W) = ™. (38) 

where v(X) is the usual density of states, on the real energy axis, for the closed chain. 
An important feature of this result is that the large-J asymptotics does not contain R. 
Therefore, by plugging in (38) into (35) we immediately obtain the simple result for the 
small- F asymptotic of the average DOR: 

<KX , Y)> ^W. m 

Note that the factors containing 1 — ^ = D 2 had cancelled, so that X enters only via 
the density of states and the localization length. This cancellation is the consequence of 
the obvious scaling property of the general equation (35), namely: the dependence on 
D drops out if one transforms to a new variable, Y — > Y/D, and to a new function, 
p — > Dp. Thus, the asymptotic (1/Y)- behavior is universal, in the sense that it holds 
for any degree of disorder and for any — 2 < X < 2. 

It remains to prove the result stated in (38). It is convenient to define the function 

oo 

W(R, J) = J dJ'P(R, J') (40) 
j 

and to study its large-J asymptotics. An integral equation for W(R, J) is derived by 
integrating (34) over the second argument of P. While performing this integration one 
should remember that the function P(R, J') is identically zero for J' < R 2 . It follows 
from this property that for J < R 2 

oo 

W(R, J)= J dJ'P(R, J') = P{R). (J < R 2 ). (41) 

o 

Combining both cases, i.e., J < R 2 and J > R 2 , one obtains the following integral 
equation: 

W(R,J) = e(l-^)P(R) + ^e(^-l)J dR' q (x-±-R')W(R',jL-l). (42) 
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We are interested in the large- J behavior of W(R, J). It is useful to define the Laplace 
transform 



oo 



W(R, s) = J dJe~ sJ W(R, J) (43) 
o 

and to study its small-s behavior. The integral equation for the transform is obtained 
directly from (42): 

W(R, s) = W (R, s) + e- sR2 f dR'q(x - 4 - R')W(R', sR 2 ), (44) 

J R 



where 

W (R,s) = -a~e- sR2 )P(R). (45) 
s 

Anticipating that W(R, J) is proportional to 1/J, we look for the solution (in the small 
s limit) of (44) in the form 

W(R,s) = A(R)\ns + B(R) + --- (46) 
This Ansatz satisfies Eq.(44) if 

A(R) = f dR'q{x - 4 - R')A(R') (47) 
J R 

and 

B(R) = B (R) + f dR'q{x - 4 - R')B(R'), (48) 
J R 

with 

B (R) = R 2 P(R) + A(R) In R 2 . (49) 

Since the site energy distribution, q(e), is normalized to 1, it immediately follows from 
(47) that A(R) is a constant. We denote this constant by —a and determine it as follows. 
Define a function $(i?') which is the non-trivial solution of the integral equation 

f dR${R)q{x- 4- R')- (50) 
J R 

Note that the kernel q(x — — R') is not symmetric with respect to R and R' and that 
equation (50), unlike (47), is not solved by a constant. The solution of (50) is: 

•<*) = (51) 

We now go back to (48), multiply it by $(-R), integrate over R and, with the help of (50), 
obtain: 

J dR&(R)B (R) = 0, (52) 
10 



which, using (49) and the aforementioned result A(R) = const = —a, yields the value of 



a: 



J dRP(R)P(^) 
a ~ J ^P(l)ln^> ( 53 ) 

where the explicit form of (see (51)) has been used at the last step. The integral 

in the numerator is known in localization theory [12] and is equal to the density of states 
z/(X). The integral in the denominator is equal to (2/£(X))(see (37)). Thus, in the small- 
s limit, W(R, s) = —a Ins, with a = v{X)^{X)/2. As a consequence, in the large- J limit, 
W(R, J) = a/J and P(R,J) = a/J 2 which completes the proof of (38). One can also 
show that (48) has a unique solution B(R) and the behavior of W(R,s) suggests that 
W(R, J) = a/J + R 2 B(R)/J 2 + ... and p(R, J) = a/J 2 + 2B(R)R 2 / J 3 + ... 



6 Conclusions 

We studied the average density of resonances (DOR) for a semi-infinite disordered chain, 
coupled to the outside world by a (semi-infinite) perfect lead. The main result of this 
work is the expression (18) for the average DOR, supplemented by the three integral 
equations, (22-24), for the three functions, h, g and /. This set of equations provide the 
general framework for calculating the average DOR, for an arbitrary disorder and coupling 
strength. We applied these general equations to the case of weak coupling and derived 
an asymptotically exact expression for the average DOR, in the limit of small resonance 
width (Eq.(39)). This expression is universal, in the sense that it holds for any degree of 
disorder and everywhere in the (unperturbed) energy band. 

It is worthwhile to emphasize the essential difference between the average DOR, as 
defined in this paper, and the probability distribution of the resonance widths, often used 
in the literature (e.g. [5]). For a finite size chain, of N sites, the two quantities differ only by 
a factor N. However, when iV approaches oo, the probability distribution shrinks towards 
a delta- function, while the average DOR approaches a perfectly well defined limit. Since 
we have taken a semi-infinite chain from the start, it was essential for us to work with the 
average DOR, rather than with the probability distribution of the resonant widths. This 
enabled us to obtain the 1/F-behavior in the N — > oo limit, in contrast to the statement 
of [5] that in that limit the 1 / V-behavior has no region of applicability. 
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